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ABSTRACT 

We report the discovery of the second accreting millisecond X-ray pulsar (AMXP) in the globular 
cluster NGC 6440. Pulsations with a frequency of 205.89 Hz were detected with the Rossi X-Ray 
Timing Explorer on August 30th, October 1st and October 28th, 2009, during the decays of < 4 
day outbursts of a newly X-ray transient source in NGC 6440. By studying the Doppler shift of the 
pulsation frequency, we find that the system is an ultra-compact binary with an orbital period of 57.3 
minutes and a projected semi-major axis of 6.22 light-milliseconds. Based on the mass function, we 
estimate a lower limit to the mass of the companion to be 0.0067 M Q (assuming a 1.4 M Q neutron 
star). This new pulsar shows the shortest outburst recurrence time among AMXPs (~ 1 month). If 
this behavior does not cease, this AMXP has the potential to be one of the best sources in which to 
study how the binary system and the neutron star spin evolve. Furthermore, the characteristics of 
this new source indicate that there might exist a population of AMXPs undergoing weak outbursts 
which are undetected by current all-sky X-ray monitors. NGC 6440 is the only globular cluster to 
host two known AMXPs, while no AMXPs have been detected in any other globular cluster. 
Subject headings: 



1. INTRODUCTION 

Accreting millisecond X-ray pulsars (AMXPs) are 
rapidly spinning neutron stars accreting matter from a 
low-mass stellar companion. These systems are thought 
to be the progenitors of the millisec ond radio pulsars 
(|Afpar et all 119821; iBacker et all 119821 ) . Therefore, their 
study allows us to better understand how neutron stars 
in low-mass X-ray binaries (LMXB) evolve into radio 
millisecond pulsars. Furthermore, AMXPs are perfect 
laboratories to study the interaction between the neu- 
tron star magne tosphere and the accretion disk (see, e.g., 
lPoutanenll2006l and references therein). 

A total of ten AMXPs among the ~ 100 neutron 
star LMXBs were known by August 2009. Seven of 
these AMXPs showed outbursts lasting a few days to 
months, and recurrence times of at least 2 years. 
In all seven cases, the pulsations were detected per- 
sistently during the whole outburst. The remaining 
three AMXPs showed pulsations only i ntermittently 
(iGallowav et all 120071; iGayriil et all 120071: ICasella et all 
120081: lAltamirano et all 120081: iPatruno et all 120091 ) 
bridging the gap between the small number of 
AMXPs and the large group of non-pulsating LMXBs 
(|Altamirano et al. 2008). The reason why only a small 
amount of neutron star LMXBs show pulsations is still 
under debate. One of the proposed scenarios for the 
non-pulsating systems is that the neutron-star magnetic 
field could be temporarily buried by the accreted matter 
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(ICumming et alll2001[) . If the time -averaged mass accre- 
tion rate of the accretors is relatively high, the accreted 
matter can bury the field throughout the life of the X-ray 
binary so the neutron star does not pulsate. However, if 
the average accretion rate is low enough that the mag- 
netic field can diffuse through the accreted matter faster 
than it is buried, then these systems will probably ex- 
hibit pulsations. Another of the proposed scenarios that 
explains why the majori ty of neutro n stars in LMXBs 
do not pulsate is that of lLamb et all (j2009f ): if neutron 
stars undergo long periods of accretion (at high rates), 
then their magnetic poles are naturally forced to align 
to their spin axes as they spin up. When the accre- 
tion rate decreases, the magnetic poles can move away 
from the rotation axis (due to magnetic dipole and other 
braking torques which cause neutron stars to spin down) 
and oscill ations powered by accre tion should become visi - 
ble. Both lCumming et all (|2001[ ) and lLamb etal] (120091 ) 
present clear predictions: we should find pulsations in 
many (if not all) of the systems which show low time- 
averaged mass accretion rates. 

Very recently, a Chandra X-ray observation serendipi- 
tously discovered a new X-ray transient in the globular 
cluste r NGC 6440 (IHeinke et all 2009c; Heinke fc Budacl 
2009; Heinke et al.ll2009bl identified as NGC 6440 X-2). 
RXTE and Swift follow-up observations not only showed 
that this source has a very low ti me-average accretion 
rate (< 2 • 10~ 12 M Q year -1 , see IHeinke et all [2009al 
but, as predicted by the m odels discussed above, al so mil- 
lisecond X-ray pulsations (lAltamirano et all 120 091 One 
of the most interesting properties of this new AMXP 
is its very short re currence time (~ 1 month, see also 
IHeinke et al.ll2009al) . If its outbursts continue, it will be 
possible to study the evolution of spin and orbital pa- 
rameters on timescales never before probed. 

In this Letter we report the discovery of X-ray pul- 
sations in NGC 6440 X-2. In a companion work 
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Fig. 1. — Intensity (2.0-10.0 keV cts/sec/5 PCU) vs. time as mea- 
sured by the PCA Bulge Scan monitoring (Swank & Markwardt 
120011 ') . The dashed lines mark the times when we detected X-ray 
pulsations (August 30th, October 1st and October 28th, 2009); 
not e that the bul g e scans did not sample the October 1st outburst. 
Sec Hcinkc ct al. (2009a) for more details. 



(jHeinke et al.l l2009ah . we report on a detailed descrip- 
tion of the evolution of the source's outbursts based on 
our multiwavelength campaign on this source. Our re- 
sults demonstrate that there may exist a population of 
AMXPs undergoing low-luminosity outbursts which are 
undetected by current all-sky X-ray monitors. 

2. LMXBS IN THE GLOBULAR CLUSTER NGC 6440 

NGC 6440 i s a glo bular cluster at 8.5 ± 0.4 kpc 
(lOrtolani et al.l I1994D . Using Chandra images, 
IPoolev et alH 2002 N ) studied the population of the low- 
luminosity X-ray sources in NGC 6440. These authors 
found 24 X-ray sources above a limiting luminosity of 
~ 2 • 10 31 ergs s _1 (0.5-2.5 keV) inside the half-mass ra- 
dius of the cluster (all of which lie within ~2 core radii 
of the cluster center); there is strong evide nce that 8 of 
these systems are LMXB s in quiescence (jHeinke et al.l 
2003). IPoolev et all (|2002h also found excess emission in 
and around the core of NGC 6440, suggesting unresolved 
point sources. 

Bright X-ray outbursts from a LM XB were ob- 
serve d in 1971, 1998 2 001 a nd 2005 dMarkert et all 
1971 lm^ Zand et all Il99l iVerbunt et al l EOOOT 
in't Zand et al. 200li: iMarkwardt fc Swankll2005l) . From 
X-ray and optical observations, lin't Zand et all ()2001l ) 
concluded that the 1998 and 2001 outbursts were from 
the sa me object, which t h ey de signated SAX J1748.9- 
2021. lAltamirano et all (|2008[ ) detected intermittent 
X-ray pulsations at ~ 442 Hz in the 2001 and 2005 
outbursts (see also iGavriil et al.l [2007h . Whether the 
1971 outburst was from SAX J1748.9-2021 is uncertain. 

Between MJD 54875 and 55075 at least 5 short- 
duration (~days) outbur sts have been dete cted in 
NGC 6440. As discussed in lHeinke etHI (|2009al) . Chan- 
dra, RXTE and Swift data indicates that all outbursts 
were from the same source, denoted NGC 6440 X-2. 

3. OBSERVATIONS, DATA ANALYSIS AND RESULTS 

We used data from the RXTE Proportional 
Co unter Array (PCA , for instrument information 
see iJahoda et al.l 120061) . Between August 1998 and 
November 4th, 2009, there were 49 pointed observations 



TABLE 1 
Timing parameters for NGC 6440 

Parameter Value 

Orbital period, P OT .f, (seconds) 3438(33) 

Projected semi major axis, a x sini (lt-ms) . . 6.22(7) 
Time of ascending node, T asc (MJD/TDB) 55073.0344(6) 

Eccentricity, c (95% c.J.) < 0.07 

Spin frequency v (Hz) 205.89215(2) 

Pulsar mass function, f x (xlO~ 7 Af ) 1.6(1) 

Minimum companion mass 1 , M c (Mq) 0.0067 

All errors are at Ax 2 — 1. 

1 : The minimum companion mass is estimated assuming a neutron star 
mass of 1.4 Mq. 

of the globular cluster NGC 6440 each covering 1 to 5 
consecutive 90-min satellite orbits. Usually, an orbit 
contains between 1 and 5 ksec of useful data separated by 
1-4 ksec data gaps due to Earth occultations and South 
Atlantic Anomaly passages. The first 27 observations 
sample 3 out bursts of the Intermitt ent AMXP SAX 
J1748.9-2021 (|Altarnirano et al.l l2008h . The remaining 
22 were triggered based on results of the PCA Bulge 
Scan program or ToO Swift observations, and sample 
the 5 outbursts of NGC 6440 X-2 (see also lHeinke et all 
I2009al) . 

We performed a timing analysis using the 
high-time resolution data collected in the Event 
(E_125us_64M_0_ls), Good Xenon and Single-Bit modes. 
Fourier power spectra were constructed for each ob- 
servation, using data segments of 128, 256 and 512 
seconds and with a Nyquist frequency of 4096 Hz. No 
background or dead-time corrections were made prior 
to the calculation of the power spectra, but all reported 
rms amplitudes are background corrected; deadtime 
corrections are negligible. 



••• 



iA 



r \r v 

0.2 0.4 06 0B 1 1.2 1.4 16 18 2 



205.75 205.8 205.85 205.9 205.95 

FREQUENCY (Hz) 



Fig. 2. — Leahy normalized (Leahy ct al. 1983) power spectrum 
of August 30th observation (~ 3 ksec of data). The signal is dis- 
tributed on at least three independent frequency bins of size 1/512 
Hz. Inset: Pulse profile in the 2—14.8 keV range (two cycles are 
plotted for clarity). The fractional rms amplitudes of the fun- 
damental and first overtone are 7.5 ± 0.6% and 2.3 ± 0.6% rms, 
respectively. Count rates (in counts/bin) are corrected for back- 
ground. 



3.1. Light curves and power spectra 
Figure Q] we plot the intensity (2.0-10.0 



keV 



In Figure yj we 
cts/sec/PCA) of the globular cl uster NGC 6440 as 
meas ured with PCA Bulge Scans (jSwank fc Markwardt! 
1200 ll) . The data suggest that between MJD 54875 and 
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55075 the globular cluster has brightened for short pe- 
riods (~days) on ~ 5 occasions (see also iHeinke et al.l 
l2009al) . However, we note that some of these points 
might be spurious given that the intensity measured is 
systematically uncertain due to the contributions from 
diffuse galactic emission and other nearby sources in the 
PCA field of view. In none of the pointed nor the Bulge 
Scan observations did we detect thermonuclear X-ray 
bursts. 

As can be seen in Figure [TJ on August 30th, 2009 
(MJD 55072) the PCA bulge scan detected the highest 
intensity in the direction of NGC 6440 over a period of 
more than 250 days. A ToO RXTE observation (Ob- 
sID: 94044-04-02-01) found NGC 6440 at a luminosity 
of L2-iokcV ~ 1-7 • 10 36 ergs s " 1 (using a distance of 
8.5 kpc, see lOrtolani et al] |l994). The energy spectrum 
is well fitted with an absorbed power law with index of 
1.83 ± 0.3 (x 2 /dof = 0.84 for 47 degrees of freedom; the 
interstellar absorption Nh was fixed to a value of 5.9- 10 21 
cm , consistent with the clu ster value and the a bsorp- 
tion measured by Chandra, see lHeinke et a l. 2009a). The 
power spectrum of this observat ion shows two featu res 
with characteristic frequencies (jBelloni et al.1 l2002t ) of 
7.4 ± 0.8 and 1.13 ± 0.06 Hz, quality factor Q of 1.3 ± 0.6 
and 0.71 ± 0.07, and rms amplitudes of 28 ± 3% and 
42 ± 2%, respectively. We searched for kHz QPOs, but 
found none. The average (background-subtracted) count 
rate during this observation was ~ 22 cts/sec and only 
one proportional counter unit (PCU) was active (2-60 
keV). Additional RXTE observations were performed on 
the 1st and 2nd of September; fluxes were consistent with 
background emission. These measurements are consis- 
tent with the fluxes measur ed by Swift between the 1st 
and the 4th of September (jHeinke et al.ll2009al) . If we 
take into account that a PCA bulge scan observation did 
not detect significant flux in the direction of NGC 6440 
on August 28th, we can estimate an outburst duration 
of no more than ~ 3 days (above ~ 10 35 erg s _1 ). 

On July 30th, October 1st and October 28th, 2009, 
RXTE sampled thr ee other outbursts of NGC 6440 
(jHeinke et al.ll2009al ) for - 1.5 - 2.5 ksec. We did not 
find any significant QPOs nor broad band noise in their 
power spectra, probably due to poor statistics (note 
that the total count rates measured are dominated by 
background counts). These observations ranged from 
^2-iokcV = 1 x 10 35 to 5.2 x 10 35 ergs s _1 , thus provid- 
ing poorer statistics than on August 30th. We found no 
evidence for dips or eclipses in any of the RXTE pointed 
observations of NGC 6440 X-2. 

3.2. Pulsations 

Adopting a source position a — 17 48 m 52 s .75, 6 = 
— 20°2 1 24 .0 (from Chandra images, see IHeinke et al.l 
l2009al) . we converted the photon arrival times to the So- 
lar System barycenter with the FTOOL faxbary, which 
uses the JPL DE-405 ephemeris along with the space- 
craft ephemeris and fine clock corre ctions to provide a n 
absolute timing accuracy of 3.4 fj,s (pahoda et al.| [2006'). 

We found a strong signal at a frequency of ~ 205.89 Hz 
in the August 30th RXTE observation. In Figure [2] we 
show the Leahy normalized power spectrum. We calcu- 
lated a dynamical power spectrum with a 512 sec sliding 
window and found that the pulse frequency is modulated 
in a way that is typical from Doppler shifts due to the 
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Fig. 3. — Upper panel: Dynamical power spectrum of August 
30th observation with a 512 sec sliding window. The frequency is 
modulated as expected from the Doppler shifts due to the orbital 
motion of a pulsar in a binary system. The data sample a little less 
than a full orbital cycle. The filled line represents our best orbital 
solution. Lower panel: Phase residuals as estimated from 300 sec 
data segments. All errors are at A\ 2 = 1. 

orbital motion of a pulsar in a binary system (see Fig- 
ure [3]) ; the data sample a little less than a full orbital 
cycle. 

To improve the signal to noise of the pulsations, we 
first modeled the pulse drifts with a sinusoid (assuming 
zero eccentricity). We found that the frequency drift is 
consistent with a system with an orbital period of 57 min- 
utes and a projected semi-major axis of 6.2 lt-ms. With 
this provisional solution, we folded the 3 ks of lightcurve 
into 10 pulse profiles of ~ 300 s each and fitted them with 
a constant plus two sinusoids representing the pulse fre- 
quency and its first overtone. We then phase-connected 
the pulse phases by fitting a constant pulse frequency 
plus a circular Kepleri an orbital model. Th e procedure 
is described in detail in lPatruno et al.l (|2009ft . In Table Q] 
we report the best fit. In the inner panel of Figure [2] 
we show the folded pulse profile in the 2-14.8 keV range 
(i.e. channels 0-35) after correcting for the effects of 
the orbital modulation. The fractional rms amplitudes 
of the fundamental and first overtone are 7.5 ± 0.6 and 
2.3 ± 0.6%, respectively (2.0-14.8 keV). 

We applied the above technique to 
bands: 2-5.7 keV, 6.1-14.8 keV and 
(We note that we could only choose 
as the data were split in two Single-bit modes - 
SB.125us.0_13.ls & SB.125us_14_35.1s - and one Event 
mode - E_62us_32M_36_ls.) The pulsed fractional rms 
amplitudes of the fundamental were (7.5 ± 1.0)% and 
(8.2 ± 1.0)% in the 2-5.7 keV and 6.1-14.8 energy band. 
The first overtone was significantly detected in the latter 
band, with an amplitude of 3.3 ± 0.8% rms. Due to the 
low quality of the data at high energies, only upper lim- 
its of 21% rms (at 95% confidence level) can be put on 
the fundamental in the 15.2-60 keV band. (NGC 6440 
X-2 was more than 100 time s brighter than the other 24 
X-ray sources in NGC 6440 (jPoolev et al.ll2002ft . so their 
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contributions can be ignored. However, Galactic Ridge 
emission does provide additional background at this lo- 
cation, not included in standard RXTE background esti- 
mates, and therefore the rms amplitudes we quote above 
might be slight underestimates). 

We also searched for pulsations in the observations on 
July 30th, October 1st and October 28th, 2009, by cor- 
recting for the orbital motion of the system and folding 
the data around the best spin period determination. No 
pulsations were found on July 30th, with 95% c.l. upper 
limits of 8% rms on the pulse fractional amplitude. On 
October 1st and 28th we found 3.4cr and 8.9cr detections, 
respectively, of pulsations with characteristics consistent 
with those we found on August 30th. 

RXTE observed NGC 6440 during the three out- 
bursts of the AMXP SAX J1748.9-2021 (see, e.g., 
lAltamirano et all 120081: IPatruno et al.ll2009t ). Although 
improbable, it is possible that NGC 6440 X-2 was also ac- 
tive during these periods. To further investigate this, we 
searched all previous RXTE observations of NGC 6440 
for the presence of pulsations at ~ 205 — 206 Hz. We 
barycentered all data in the same manner as above and 
divided the data in 3 groups, corresponding to the 1998, 
2001 and 2005 outbursts. For each group we performed 
Fourier transforms of all the data and calculated a single 
averaged power spectra. We found no significant pulsa- 
tions in the ~ 205 — 206 Hz range. We also searched for 
evidence of inte rmittent pulsations such as those found in 
other AMXPs ( Gallow ay et"alll2007t IGavriil et all 120071 : 
lAltamirano et all 120081: ICasella et al1l2008f ). but found 
none. In the above searches, we did not correct for the 
orbital motion of the binary; although these corrections 
would improve our sensitivity, such analysis is beyond 
the scope of this Letter. 

Finally, we note that it was not possible to search for 
pulsations in the Swift observations. These data-sets 
were obtained in photon-counting mode with a time res- 
olution of 2.5 seconds. 

4. DISCUSSION 

We have discovered the second accreting millisecond 
X-ray pulsar in the globular cluster NGC 6440. Pulsa- 
tions were detected in three RXTE pointed observations 
performed during the decay phase of thr ee outbursts, 
each o f which lasted less than ~ 3 days (|Heinke et al.l 
l2009al) . The frequency drifts we observe are consistent 
with an orbital period of 57.3 minutes and a semi-major 
axis (a x sini) of 6.22 light-msec, which shows that this 
new AMXP is in an ultra-compact binary system. 

The characteristics of this new pulsar are similar to 
those of other AMXPs, particularly SWIFT J 1756.9- 
2508 (jKrimm et al.l 120071: IPatruno et all I2010D . which 
showed pulsations at a frequency of 182 Hz, an orbital pe- 
riod of 54.7 minutes, a projected semi-major axis of 5.94 
lt-ms and a min imum companion m ass of 0.0067 Mq. As 
pointed out bv lKrimm et al.l 120071 such a low mass for 
a companion in this type of system is probably inconsis- 
tent with brown dwarf models while white dwarf models 
suggest that the companion is probably a He-dominated 
donor. 

The reason why only twelve 5 NS LMXBs exhibit coher- 
ent pulsations out of about 100 systems is still unclear. 



There is a strong tendency for the detected AMXPs to 
have rather low time- averaged accretion r ates (over many 
outbursts, see, e.g., iChakrabartvl 120051 iGallowavl 12006: 
IWiinandsl [20081 ) . The time -averaged accretion rate of 
NGC 6440 X-2 is at most similar to that of the AMXP 
SAX J1808. 4- 3658, but could easily be much lower (see 
discussion by iHeinke et al.l l2009af) . This is consistent 
with the general tend ency and furthermore, with the the- 
oretic al predictions (|Cumming et al.l 120011 : lLamb et al.1 
2009) that a significant fraction (if not all) of the NS 
LMXB systems with low time-averaged accr etion rates 
harbo r accreting millisecond X-ray pulsars (W iinandsl 
l2008h . 

The outburst durations (~ 3 days) of the newly 
discovered AMXP are short compared with those of 
the other AMXPs, although not uniquely so. For 
example, short-duration outbursts of the AMXP 
XTE J 175 1 -305 have been dete c ted in 2005 2007 
and 2009 dGrebenev et all 120051: iSwank et al.l " 
Markwardt fc Swankl 120071 IMarkwardt et al.l 



2005; 



20071 



Linares et al.l l2007t IMarkwardt et al.l l2009bl ). Our dis 



covery of NGC 6440 X-2 indicates that there might exist 
a population of AMXPs undergoing weak outbursts. 
Unfortunately, short low-luminosity outbursts cannot 
be detected with the current all-sky X-ray monitors 
(e.g., ASM aboard RXTE). More sensitive instruments 
like the "Monitor of All-sky X-ray Image" (MAXI, 
which should detect X-ray - 2-30 keV - sources of 
about 20 mCrab from 90 minutes observations, and 
sources of about 4.5 m Crab after a day observation, see 
iMatsuoka et al.|[2009D will probably allow us to discover 
many more poten tial low-luminosity A MXPs. 

As discussed in IHeinke etHl (j2009al ). NGC 6440 X-2 
has been in outburst at least five times in the last 250 
days. If true, this is the AMXP with the shortest re- 
currence time. If NGC 6440 X-2 continues undergoing 
this type of short outbursts and RXTE is able to ob- 
serve them, then NGC 6440 X-2 has the potential to be 
the best millisecond X-ray pulsar in which to study the 
evolution of the binary as well as that of the neutron star 
spin (by means of coherent timing analysis). 

The globular cluster NGC 6440 is known for having 
at least 24 faint X-ray sources to a limiti ng luminosity 
of ~ 2 • 10 31 ergs s" 1 (0.5-2.5 keV, see iPoolev et ail 
120021 ). NGC 6440 is also known for hosting at least 
6 millise cond radio pulsar s, of which 3 are in binary 
systems (|Freire et al.l 12008). None of the known radio 
pulsars match the position (nor the characteristics) of 
NGC 6440 X-2. Although millisecond radio pulsars in 
binary systems have been found in o ther globula r clus- 
ters (see, e.g., iRansom etaLl 120051: iFreire et all l2008l 
and references therein), it is an intriguing question why 
NGC 6440 is the only globular cluster kn own today to 
host (two ) AMXPs (SAX J1748 9-2021 - IGavriil et all 
[20071 and lAltamirano et al.ll2008l - and NGC 6440 X-2 
- this paper), while no AMXPs have been detected in 
any other globular cluster. This might be only a se- 
lection effect, as current all-sky monitors cannot detect 
short low-luminosity outbursts. Monitoring X-ray obser- 
vations of different globular clusters would be needed to 
further investigate this. 
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5 See also the rece nt discovery of a 245 Hz AMXP by 
IMarkwardt et al.1 ((20093). 
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